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Abstract. A general scheme is developed to deal with ID lattice systems that could 
be topologically complicated. It is aimed to give a complete study of two coupled 
normal metal rings. Our method starts with an investigation of the local expressions 
of the eigenfunctions. By connecting different parts of the system, all the eigenvalues 
and eigenfunctions can be obtained. It is found that there is a possibility for the 
existence of localized states, which is beyond previous expectations. 
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1. Introduction 

The tight binding model is widely used to study electron behaviours in crystals. While 
it is originally suited for real materials that are of 3D, people are tempted to study 
the ID case for simplicity. Usually only bulk properties are considered, so the Born- 
von Karman periodic boundary condition is chosen for convenience. The dispersion 
e = — 2t cos/c is well-known for this ID periodic tight-binding model [T]. Although 
open boundary conditions are most encountered in reality, the argument often used is 
that the physics in the bulk should not depend on any boundary conditions within the 
thermodynamic limit. 

As technology develops, ever smaller systems can be manufactured, and realistic 
ID systems like quantum wires are realized. It happens that certain systems become 
sensitive to their boundary conditions due to small sizes. Additionally, the ID models 
that theorists could only play with previously come to be relevant. In this background, 
some combined ID systems have been proposed, for example two coupled rings 0 0. 
To deal with models with complicated boundary conditions, a systematic method is in 
need. 

In previous works on the two coupled rings, several different methods have been 
used. In [2], the method of Dirac constraints [I] was used. In j3], a simple wave function 
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ansartz method was used . Since the focus of these papers was on persistent current 0® 
[7], the solutions of the coupled system were discussed only briefly and not completely. 
The present work is dedicated to outline a general scheme to solve complicated ID 
systems. By “solve”, it is meant to obtain all the energy levels and the corresponding 
eigenfunctions. The two coupled rings problem will be solved completely under this 
scheme. 

Instead of the continuous model which is used in 0 0 . a lattice model is applied 
in this work. The physics has no essential difference between these two models, but the 
latter one has two advantages. First, the lattice model follows exactly the spirit of tight 
binding, such that when the system becomes small to the extent of several nanometers, 
the discreteness of atoms may become important. Second, the results could be checked 
easily by diagonalizing Hamiltonian matrices numerically. 

The usual way of doing quantum mechanics in ID is to divide the system into several 
homogeneous parts, obtain the general form of the wave function for each part first, and 
then match them at the boundaries [8]. We will employ this standard procedure to study 
two coupled rings. Before we jump into this complicated problem, it is pedagogical and 
beneficial to start from simple situations with simple boundary conditions. Thus we 
would like to give a systematic review of ID tight binding models, and see how different 
boundary conditions could be treated. 

We will begin with the semi-infinite wire in Section 2. Here the single open end of 
the wire is the natural starting point for us to obtain the eigenfunction by iteratively 
applying the Schroedinger equation. We will use a transfer matrix method to do this 
iteration. The general form of the eigenfunction is determined by calculating the power 
of the transfer matrix. Since there is only one boundary condition for the semi-infinite 
wire, the eigenfunction is obtained readily. After this model is fully understood, we will 
consider more complicated boundary conditions. 

In Section 3 we study a finite open wire with two open ends. We use the same 
scheme by starting from one end and using the transfer matrix to obtain a general form 
of the eigenfunction, but now the eigenfunction should also fulfil the boundary condition 
on the other open end. This gives us an additional equation which gives a constraint 
to the eigenvalues and so determines a discrete spectrum. In contrast, the semi-infinite 
wire has no such equation and it has a continuous spectrum. 

In Section 4 we come to study a single closed ring where the boundary condition is 
periodic. The result for this case is well known, since a Fourier transform can be applied 
to solve it very quickly. While now we don’t want to use the Fourier transform method, 
we want to perform the same procedure as we have done for the finite open wire. This 
offers us an opportunity to test our scheme by comparing our results to the familiar 
ones. The boundary condition now is distinct from the finite open wire, so we have to 
obtain a different equation that determines the eigenvalues. An interesting point for 
the single closed ring is that we can add magnetic flux through it, which leads to the 
famous Aharonov-Bohm effect [5] and persistent currents |7]. For this situation, a gauge 
transform can be used to give the general result. 
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In section 5 we will begin to study two coupled rings. Here we assume the rings 
are connected by allowing electrons to hop from a site on one ring to a site on the other 
ring. It is different from the model considered in 0 0. in which there is a common 
site for both rings. Actually our current scheme is capable to deal with both situations. 
But since our goal is only to outline a general scheme and it is not aimed to study all 
the possible models, we do not consider the common-site model in this work. To solve 
the problem of the two coupled rings, first we get the local form of the eigenfunction for 
each single ring, just as what we have done for the single closed ring. Then we use the 
boundary conditions to connect the local wave functions, where the connection equations 
are unambiguously written down by observing the Hamiltonian matrix directly. Similar 
to the finite open wire and the single closed ring, we obtain an equation that determines 
the complete spectrum of the system. By studying this equation carefully, we find there 
may exist localized states. The possibility of this fact is ignored in previous works [2j 
0- 


2. Semi-Infinite Wire 


We use the single band tight-binding model throughout this paper. It means on each 
lattice site there is one and only one state for a single electron. An electron can hop from 
one site to its nearest neighbours. Electrons are assumed spinless and non-interacting. 
The many particle states are composed by occupying single electron levels one by one. 

Our starting point is a semi-infinite wire, of which the configuration is shown in 
Figure [lj The most notable feature in this figure is that there is a single open end for 
the semi-infinite wire. The Hamiltonian under our assumption can be written in the 
formalism of second quantization as 

OO 

H = —t ^ a\ +1 ai + h.c. (1) 

i=0 

where t is the hopping constant, a\ and a, are creation and annihilation operators for 
electrons on site i. 

Set A = (ao,ai, •••) T , and rewrite the Hamiltonian in a matrix form 

H = -tA^hA ( 2 ) 


where 


h 


(0 1 \ 
1 0 1 
1 0 1 
1 0 

V ••• / 


(3) 


is equivalent to the Hamiltonian matrix in the formalism of first quantization. Assume 
0 = (xo, Xi, ...) T to be an eigenvector of h with eigenvalue A, then [x 0 ao + X\a\ + x 2 a\ + 
■ ■ -]|fi) (|Q) represents the vacuum state) is a single particle eigenstate of H with energy 
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Figure 1. Lattice model for the semi-infinite open wire. The solid spots on the line 
represent lattice sites. They are linked by allowing electrons to hop between nearest 
neighbour sites. 

level e = —At. We would like to diagonalize h, that is to say, we want to find all the 
eigenvalues of h and the corresponding eigenvectors. In this way we can obtain all the 
single particle eigenstates of H. 

The eigenvalue equation of h reads 


(0 1 \ 


f Xo \ 


{ x ° \ 

1 0 1 


X\ 


Xl 

1 0 1 


x 2 

= X 

X 2 

1 0 


X 3 


X 3 

V •••) 


V : ) 


V : / 


It is also the time independent Schroedinger equation in the formalism of first 
quantization. Write equation Q explicitly 

X\ = Xx 0 
xo + x 2 = Aaq 
xi + x 3 = Xx 2 
x 2 + x 4 = Xx 3 


Now if we set xq = 1, then x\ = A. For n > 1, x n — Xx n -\ — x n - 2 . By combining an 
obvious identity x n _i = i„_i, we have 


where 




-1 

0 


(5) 

( 6 ) 


is the so-called transfer matrix. 

L is a constant matrix, independent of n. We can iterate equation ([5]) and express 
(. x n ,x n -i) T in terms of (xi,a:o) T as 


n 1 


= L n 


Xi 

Xq 


(7) 


Now if we know L n we can obtain the entire wave function. The calculation of L n 
depends on the value of A. We summarize the results here and put its derivation in 
Appendix A. 










5 


i) if A G (—2, 2), set A = 2 cos 9 , where 9 G (0,7r), then 


L n = 


sin(n+l)0 
sin 9 
sin nO 
sin 6 


_sin nO 

sin 9 

sin(n— 1)9 
sin 9 


ii) if A = 2, then 


L n = 


n + 1 — n 


n — (n — 1) 

which corresponds to 9 —>■ 0 in (|8]) . 
if A = —2, then 

n + 1 n 
—n —{n — 1) 

which corresponds to 9 —y tt in 


L n = (-l) r 


hi) if A > 2, set A = 2coshfc, where k > 0, then 

sinh nk 


L n = 


sinh(n+l)fc 
sinh k 
sinh nk 
sinh k 


sinh k 

sinh(n— l)k 
sinh k 


which corresponds to (|8j) if we set 9 — ik 

if A < —2, set A = —2 cosh k, where k > 0, then 

t sinh(n+l)fc sinh nk 

sinh k sinh k 

sinh nk sinh(n-l)fc 

sinh k sinh k 

which corresponds to (pi) if we set 6 = 7r + ik 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


Once we know L n , we can plug L n into ([7]) and use the initial value (xo, x\) = (1, A) 
to write down the eigenfunction for each case as above: 


i) if A = 2 cos 9 G (—2,2) , 9 G (0,7r) 

, . „ sin 36* SU146 1 

^ = (1,2 cos 9, 


sin 9 ’ sin 9 ’ 

Note that £4 (cos 6*) = is the second kind of Chebyshev polynomial [9J. 


ii) if A = 2, 


if A = -2, 


"0 — (1) 2, 3,4, • ■ -) : 


0 = (l,-2,3,-4,---) n 


(13) 


(14) 

(15) 


iii) if A > 2, A = 2 cosh k, k > 0 


, , , , sinh3fc sinh 4k 

4=1,2 cosh k, • • • T 

sinh k sinh k 


( 16 ) 
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if A < - 2 , A 


= —2 cosh k, k > 0 


0 = ( 1 , —2 cosh k 


sinh 3 k 
’ sinh k ’ 


sinh 4 k 
sinh A; ’ 



(17) 


In quantum mechanics wave functions are required to be normalizable, so they 
cannot blow up at infinity. Solutions ii) and iii) therefore should be abandoned. Thus 
the spectrum for a semi-infinite wire can only be in the range of (— 2 , 2 ). Since we have 
no other constraint to the eigenvalue now, any value in (— 2 , 2 ) belongs to the spectrum, 
so the spectrum of the semi-infinite wire is continuous. 

In i), if we set xo = sin#, the eigenfunction becomes 

0 = (sin 9, sin 29, sin 39, sin 49, ■ ■ •) T (18) 

This form of the eigenfunction looks like a standing wave with only one fixed end. It 
can be considered as the superposition of the initial wave and the reflected wave when 
an electron travels from infinity to the boundary and then reflects back. 

Now we have finished the discussion of the semi-infinite wire. The key to the 
solution in this section is that we can express the entire wave function in terms of 
(x 0 ,Xi) T through ([7]). We note that the simple geometry of the semi-infinite wire gives 
us two great advantages. First, the single end of the wire supplies a natural starting point 
to begin the iteration. (Although in principle we can start from any two adjacent sites 
( x n , x n+1 ) T and derive the entire wave function, it makes calculation more complicated.) 
Second, once we have solved for the eigenfunction, we need only further to consider the 
requirement of normalization, since there are no other boundary conditions. 

In the following we want to vary the boundary conditions little by little, and see 
how our scheme can be adjusted readily for new situations. 


3. Finite Open Wire 

In this section we discuss the finite open wire, the configuration of which is shown in 
Figure [ 2 } Compared to the semi-infinite wire, it has two open ends and only a finite 
number of sites. Assume the number of sites is N, and denote the wave function as 

0 = (Xi, ...,X N ) T . 

We would like to follow the same steps as before and start from the open end, since 
there is no difference in the local region of the end as compared to the semi-infinite wire. 
Let’s begin with x± and X 2 and use the transfer matrix to get the value of x%, X 4 , • • 
Xjv- Now the difference from the semi-infinite wire is that we cannot go any further, 
because no site (N + 1 ) exists. We have to stop here and require xn -i = Ax at, which is 
the boundary condition at the other end of the wire. 

It’s appropriate to use a little trick at this point. Imagine that we add site 0 
at the left end and site (N + 1) at the right end, and require that Xo = Xat+i = 0. 
Then the boundary conditions x 2 = Axi and xat_i = Axjv could be tailored as 
x 2 = Axi — xq and xjv+i = Ax^v — %n-i, which fits the general form of the iteration 
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equation x n+ i = Xx n — x n -\. The advantage of this trick is that we have two much 
simpler boundary conditions x 0 = Xjv+i = 0 now. 


0 


2 3 


N-1 


N N+1 


Figure 2. Lattice model for the finite open wire. The solid spots in the middle of 
the line represent real sites and the star-like spots at the two ends represent imaginary 
sites. 

We start from xo( = 0) and Xi, use the transfer matrix and get a general expression 
of x n , for any n. The form of x n should be the same as that for the semi-infinite wire. 
Actually we can set X\ — 1, so X 2 — A, and then the same iterations as that of the 
semi-infinite wire follow. 

After obtaining the expression x n we need to require apv+i = 0. It’s obvious that 
there is no means for x n to increase as n increases, so we have to abandon solutions ii) 
and iii). The eigenvalue A stays in the range (—2, 2) just like the case of the semi-infinite 
wire. For the semi-infinite wire any numbers in (—2, 2) are possible spectrum, but now 
the additional equation xn+i = 0 is a constraint to A and selects some specific values in 
(- 2 , 2 ). 

Set A = 2 cos# (0 < 9 < 7r), and xo = 0, x\ — sin#. Applying the transfer matrix 
(J8|, we have x n = sin nd. From xjv+i = 0, there is sm(N + 1)# = 0. This equation has 
N distinct roots 

777/ 

Am 2 COS #m; #m y ~ £7T, Ul 1,2,..., N (19) 

They correspond to energy levels e m = —A m t which constitute the complete spectrum. 
The eigenfunction with the energy level e m is 

-0m ~ (0, sin # m , sin 2 9 m , ■ ■ ■, sin N9 m , 0) T , (20) 

up to a normalization constant. It represents a standing wave that has two fixed ends. 

Notice that the continuous version of our lattice model of the finite open wire is 
that of the particle in a box. From their solutions, the similarity is obvious. 

4. Single Closed Ring 

If we connect the two ends of the finite open wire by allowing electrons to hop between 
the two end sites, then it forms a single closed ring. This configuration, as shown in 
Figure [3j is equivalent to a ID lattice with the Born-von Karman periodic boundary 
condition. Due to the translational symmetry, it’s usually solved through Fourier 
transform. Now we want to develop our scheme to deal with the single closed ring 
and compare the results to the familiar ones. 

Assume there are N sites on the ring, with the wave function 0 = (x 0 , x l7 ..., Xat_i) t . 
Just like what we have done in previous sections, we can start from x 0 , X\ and use the 





Figure 3. Lattice model for the single closed ring. The solid spots on the circle 
represent lattice sites. 


transfer matrix to attain x n , for arbitrary n. Since now we have a periodic boundary 
condition, we expect that when n is larger than N, x n should be the same as x n -N. So 
we have xn = x 0 and x^ + i = x\. For 



thus 


( 21 ) 


X i 

x 0 


— L n 


X 1 

Xq 


( 22 ) 


Before we move on, it’s useful to unify all three different situations of A by enlarging 
the domain of 9 in i). As mentioned when we wrote down the expression of L n in section 
2, we can extend the domain of 9 into the complex plane, as shown in Figure [4j to include 
all different situations. In the following we will use this generalized interpretation of 9, 
and therefore L n could be simply expressed as ([8]). 


Now insert into (22) and we have 


sin(jV+l)e 
sin 6 
sin NO 
sin 6 


sin NO 
sin 0 
sin(iV-l)(9 


sin 6 


-1 



(23) 


If Xo and X\ are both zero, then x n = 0 for any n. Thus, in order to get a non-zero 
eigenfunction, we have to first make sure (23) has a non-zero solution. So 
sin(JV-t-l)g 


sin 6 
sin NO 
sin 0 


- 1 


_ sin NO 

sin 6 
sin(7V— 1)6 

sin 6 


= o 


(24) 


The roots of this equation represent the spectrum of the single closed ring. After fixing 
9 , the ratio of x\ and xq could be readily obtained through (23). The exact values of 
X\ and a: 05 U P to a globe phase, can be determined by normalization. The eigenfunction 
can be expressed in terms of x$ and x± through ® as 


sin n9 sin(n — 1)9 

X n = - — Q - 

sm 9 sm 9 


x 0 


(25) 
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Figure 4. The domain of 9 in complex plane under the generalized interpretation 
of 9. The thick U-shape line in the figure represents the domain, which unifies all 
different situations. 

We don’t give the detailed form of the eigenvalue and the eigenfunction for the 
single closed ring right now. We want to generalize the model by letting magnetic flux 
thread the ring. We will study this generalized model carefully. The solution of this 
generalized model will automatically include that of the ring with no magnetic flux as 
a special circumstance. 

Let’s add magnetic flux $ through the ring. (This makes the system a setup for the 
Aharonov-Bohm effect and the persistent current.) The well-known Peierls substitution 
can be used to include the effect of the magnetic flux into the tight-binding model. 
Simply speaking, the Peierls substitution varies the hopping constant t by a phase 0, 
such that t is substituted by fe'A The phase 0 is related to the vector potential A via 
0 = (q/h) f^ l+1 A ■ dr. We can choose a gauge to let A have the same magnitude along 
the ring, and so the result does not depend on the site number i in the expression of qb. 
Now <h = § A • dr is the total flux. If <f> ^ 0, then A ^ 0 and thus 0 ^ 0. 

After the introduction of magnetic flux, the Hamiltonian for the single closed ring 
can be written as 

V-l 

H = — te 1( ^al +1 di + h.c. (26) 

i =o 

where a at = a 0 is assumed for the ring’s topology. 

We could turn this into a matrix form as ([2]) and would like to diagonalize the 
first quantization Hamiltonian matrix /r^,. To emphasize the difference, we change the 
notation a little and set the eigenfunction of to be ( 2 / 0 , 2 / 1 , ■ • •, 2/at— i) T - 
Now the Schroedinger equation (the eigenvalue equation of h^) becomes 

y n -ie~ 1(p + 2 /n+ie 10 = A y n (27) 

where we set A = 2 cos 6 in the generalized interpretation of 9. This equation differs from 
our familiar one x n _i + x n+ \ = Xx n only by a gauge transform. If we set y n = e~ m ^x n 
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Figure 5. The spectrum of the single closed ring. Note: N is set to be 7 in these 
plots, (a) Solutions of (34). The horizontal axis represents N9. (b) Spectrum of the 
single closed ring. The spots on the curve e = — 2f cos# represent the roots jj-2ir ± </>. 
The arrows around show how the roots move as the magnetic flux increases, (c) The 
roots (^27 t — </>) are reflected to — (^27 t — </>) . (d) The spectrum after the adjustment. 
The arrows show the direction of the motion under the increase of </>. 


then they are the same. Thus, based on the result 
sin nO 


, we have 


Vn 


sind 


. e -i(«-l )* yi 


sinn 


sind 


— e-‘"*yo 


Next, the boundary condition could be expressed similar to (22) as 


V o = 


sin NO 
sin 6 


3 _ i(A r-i )0 _ sm(iV 1)6 iN(j) 


y i 


sin# 


y o 


y i = 


_sm(N + l)6 sin NO 


sin 6 


V i - 


sin 0 


■e-’^+'^so 


(28) 

(29) 

(30) 


Although (29) and (30) are equally important now, they are not in the two coupled 
rings system that we will consider in the next section. Due to the coupling, (30) will be 
substituted by a different equation. However (29) will still hold for the topology of the 
closed loop. From (29) we can express y\ in terms of 

e'^sin 6 + e - 1 ^sin(iV — 1)6* i (iv_ 1 ^ 


y i = 


sin NO 


yo 


Then plug (31) into (28) , we have 

g i {N-n)<j> g j nn 0 _|_ s j n (7\T _ n \Q 

y n = - 


-yo 


(31) 


(32) 


sin NO 

One notable property of this wave function is i/N-n = y* n , which means that the 
system remains the same if we reverse both the direction of the magnetic flux and 
the direction of the ring. The denominator sin NO may equal to zero, but it occurs only 
accidentally, corresponding to some degeneracy of the system. It could be lifted by a 
small perturbation of the flux. 
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In (27) if we set n — 0 there is 

Vn-i «-'* + yie" p = 2 cos 


(33) 


where the periodic boundary condition is used for the substitution of y_i to Dn-i- 

(34) 


Combine (32) with (33), after simplification we get 


cos NO = cos Ncft 

The solution of (34) is shown in Figure |5^a). Generally we have 


9 = — 27t ± 
N 


(35) 


where m belongs to the set of integers resulting in 9 G [0,7r]. There are total N such 
integers, which thus constitute the complete spectrum. 

Now if <ft varies, the spectrum 9 varies accordingly. But different 9 s vary in different 
directions as shown in Figure |5^b), depending on whether they have ’+’ or ’ in (35). 

We can enlarge the domain of 9 to [—7r, 7r] and reflect all the 9 s with ’ sign to 

the interval [— ir, 0] as shown in Figure [5j^c) 

m 

(36) 


—2vr 
N 


'-27T-0) = --27T + 
'N ’ N 


In this way we obtain a unified expression for 9 

„ 771 

6 = n 27T + ( ^ 


(37) 


where m belongs to the set of integers that let 9 G [—7r, 7 t], Now as (ft varies, all the 9s 
move in a consistent direction, as shown in Figure |5](d). 

Plug (37) in (32) and we have 

Vn = e im ^y 0 (38) 

To normalize the eigenfunction, we can set yo = ^=. 

Actually (37) and (38) can be easily obtained from a direct Fourier transform. So 
this justifies our scheme. 

We have finished the discussion of the single closed ring. It serves us a good 
preparation for the problem of two coupled rings. 


5. Two Coupled Rings 

The configuration of two coupled rings is shown in Figure |6j We assume the rings have 
equal sizes and electrons are allowed to tunnel from the left ring site yo to the right ring 
site Zq. Magnetic flux (fq is threaded into the left ring and magnetic flux $2 is threaded 
into the right ring. The Hamiltonian could be cast into 

'ft = ftl + ft-2 + ^coupling 


(39) 
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Figure 6. Lattice model for the two coupled rings. The solid spots on the circles 
represent lattice sites. The link in the middle represents the electron can hop between 
y 0 and zq. Magnetic flux <f>! is threaded into the left ring and magnetic flux $ 2 is 
threaded into the right ring. 


where 


N—X 

Hi — — d\ +1 dj + h.c.] 

i =0 
N -1 

n 2 = - + h.c.] 

3=0 

"^coupling = ho(Sg6o ft - ^o®o) 


(40) 

(41) 

(42) 


Hi and H 2 are two separate single closed ring Hamiltonians. H\ is for the left ring, 
H 2 for the right ring. H CO upiing represents the coupling between the two rings. In the 
expressions of Hi and H 2 above, it is assumed a at = ao and fejv = b 0 for the ring topology. 

We emphasize here that the model in 0 0 is different. In m 0 , y 0 and z 0 are 
combined to a single site. This site then connects to y 1 , Vn-i, %i and zn-i separately. 

Set the eigenfunction to be x\) = (yi,y2, Vn-i, Vo, z o, z i, ■■■,zn-i) T - First, we can 
express the wave function on the left ring in terms of y 0 and the wave function on the 
right ring in terms of zq separately, exactly the same as what we did for the single closed 
ring. So like (32), we could have 


g i(Af n)<tn s j n n Q _|_ e m ^i s in (N — n)9 

Vn = -^ ato - y ° 

sm N 6 


(43) 


Zn 


e i(Af n )</>2 gj n n Q _|_ e m<j >2 s j n (W 

sin NO 



(44) 


Next, let’s write down the Hamiltonian (39) in the matrix form 
H = -A^hA 


( 45 ) 
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where A = (aq, a 2 ,a N -i, a 0 , b 0 , h 

..., b N . 

_i) t , and 


/ 0 te 1 ^ 1 

te -1 ^ 1 



te 




te i 1 ^ 1 




te -1 ^ 1 0 

te 1 ^ 1 


h = 

te ^ te -1 ^ 1 

0 

^0 


Vq 

0 




te~ u 

t 


V 


te i</32 


\ 


te 1 ^ 2 

0 


te 1 ^ 2 


te~ i<t)2 


te l ^ 2 


te 1 ^ 2 

0 


(46) 


The rows for y 0 and Zo show the coupling equations 
y N - l te~ l4 ' 1 + yite’ 4 ’ 1 + Vi, v 0 = A y Q 
ZN-ite -1 ^ 2 + Zite 11 ^ + y Q y 0 = A Zq 


(47) 

(48) 


where A is the eigenvalue of h (the single level energy is e = —A under our notation). 


Set A = 2t cos 6 under the generalized interpretation of 0. Plug (43) and (44) into (47) 


and (48) to get 



2tsin6 l (coslV6 l — cosAh/q) —Vo sin N6 

— Vo sin NO 2t sin 0(cos NO — cos N(j) 2 ) 

In order to get a non-zero solution, we have 

2tsin0(cos Ah9 — cosAh/q) —Vo sin NO 

—Vo sin TV# 2t sind(cos NO — cos N(j> 2 ) 


= 0 


= 0 


So 


y2 NO 

—^ 5 — = (cos NO — cos N(pi)(cos NO — cos lV</> 2 ) 
4 1 2 sin" 0 


(49) 


(50) 


(51) 


Equation (51) is the main result of this paper. The solution of (51) determines the 


spectrum of our two coupled rings system. 

It is noted that N(f>i is proportional to the total flux that threads into the left ring. 
Actually N<pi = |^7r, where $0 is the flux quantum. N(f>i is the total phase that an 
electron could acquire due to magnetic flux $1 as it goes around one circle along the 
left ring. If this phase equals to 27 t, there should be no observable effect. Now since 


the only term involving 0 1 in (51) is cosAh/q, it’s obvious that if we substitute N(j>i by 


N<fii + 27t, equation (51) remains the same. Likewise for N<p 2 and the right ring. 


The two rings are coupled through the term "Hcoupimg- If there is no coupling, or to 
say Vq = 0 (so % coupling = 0), then we expect to have two separate rings. Now through 


(51) it’s clear that if Vq = 0, we have cos NO — cos Ah/q = 0 or cos NO — cos N(p 2 = 0. It 


means the rings decouple and the electron could only move on one of them. 
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Figure 7. Solutions of (51), the equation that determine the complete spectrum of 
the two coupled rings. Note: we set N = 8, Vo = t = 1 in this plot, (a) the situation 
when the eigenvalue A is in the range [—2, 2] and we set A = 2cos0 ; (b) the situation 
when the eigenvalue A is in the range [2, oo) and we set A = 2coshfc. 


At this level, we don’t try to make a complete analysis of the equation (51). It 
should be very complicated since (51) is a transcendental equation. Instead, we want 
to use an example to show how the solutions of (51) look like. 

Let N = 8, V 0 = t = 1 and choose two random fluxes 4>i and 0 2 - Figure [7] shows a 
sketch on how the solutions of (51) look with these parameters. 

In Figure [7](a), we show the situation when A is in the range [—2,2] (6 is real 
and in the range [0,7r]). The states here are propagating ones, for which the electron 
could move freely on both rings. As we will see, the wave functions of these states 
represent charge density waves. It may seem that all the eigenstates should belong to 
this kind, but that is not the case. There are 16 sites in our setup, so we should expect 
16 eigenstates in total. We could find 14 intersections in Figure [T^a) , which correspond 
to only 14 eigenvalues and eigenstates. Thus there must exist 2 eigenstates that are not 
propagating states, of which the eigenvalues are beyond the interval [—2,2], 

If A > 2, we could use the general interpretation of 6 and set it to be i k. Figure 
0b) shows there is a root in this range. So there is an eigenvalue in the range (2, oo). 
Similarly, there exists an eigenvalue in the range (—oo, —2), although we don’t show it 
explicitly here. 

The eigenstates of which the eigenvalues are not in [—2, 2] correspond to localized 
states (bound states). As we will see, in this situation the wave function is localized 
around the junction and decays exponentially as it goes away from the junction. We 
have never met such kind of localized states before. It appears neither on the semi¬ 
infinite wire, nor on the finite wire, nor on the single closed ring. Also, we note that 
the possibility of the existence of the localized states has never been mentioned in the 
previous works on two coupled rings. 
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Next we want to discuss the local density of states and the normalization of the 


wave functions. Now assume 0 is a solution of (51). From (43) and (44), we have 
12 


Wn 


\yo? 


*0 


From (49), there is 

W 

M 2 


sin 2 nO + sin 2 (lV — n)0 + 2 sin nO sin(A r — n)6 cos N<pi 

sin 2 NO 

sin 2 nO + sin 2 (iV — n)9 + 2smn6 sin(N — n)6 cos iV</> 2 

sin 2 NO 


(cos NO — cos N(j> 1 


(cos NO — cos N(j) 2 ) 2 


(52) 


(53) 


(54) 


Equation (52), 


and (54) show the relative magnitude of the eigenstate’s local 


density. When 0 is in the range [0,7r], it’s obvious that the local density is oscillating 
on each ring. So it represents a charge density wave and we say that the state is a 
propagating one. When 0 is out of the range [0,7r], we could use k to substitute 0 either 
by 0 = ik (when A > 2) or by 0 = n + ik (when A < —2). Then cos nO and sin nO could 
be substituted by ± cosh nk and ± sinh nk. Since cosh nk and sinh nk are basically 
proportional to e nk , it’s easy to see that the local density is largest at the junction 
region, and decays exponentially away from it. In this sense, we have a localized state. 

We can use the trigonometric identities 


iV " 1 N 

sin 2 nO = — 

^ 2 

n=0 


sin N0cos(N — 1)0 


2 sin 0 


N -1 


^ sin 2 (IV 


n=0 
N-l 


N 

2 


sin TV# cos (N + 1)0 
2 sin 0 


2 sin nO sin (N — n)0 = cot Osin NO — N cos NO 


n =0 


to sum (|52]) and (J53j) for all sites n 

N(l — cos NO cos Nfa) + cot0smN0(cosN(j)i — cos NO) 


N-l | 

E 1 ’ 

n=0 


^ \yo\ 


N ~ l |_ 12 
\ ^ \ z n\ 


sin 2 NO 

N{1 — cos NO cos N(p 2 ) + cot Osin NO (cos Nfa — cos NO) 

sin 2 NO 


(55) 


(56) 


We can combine (54), (55) and (56) to normalize the wave function and obtain the 


values of \y 0 \ and \z 0 \. The exact values of y 0 and z 0 are determined by equation (49) 


and the normalization of the wave function. Since the final expression is very long and 
gives little insight, we do not write it explicitly here. 
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6. Conclusion 

We have developed a general scheme to study ID single band tight binding models with 
complicated boundary conditions. The two coupled rings problem has been analysed 
step by step in this framework. The eigenvalues and eigenfunctions are given in a 
complete way. 

Our scheme starts from an investigation of the local property of the wave function. 
We find that in a homogeneous region, an explicit form of the wave function can be 
readily obtained in terms of its values on two nearest neighbour sites by a transfer 
matrix method. Then with the boundary conditions for some special sites expressed 
specifically, we are able obtain an equation to determine the complete spectrum. 

Note that we divide the range of the eigenvalue A mainly into two different cases. 
One is in [—2, 2], the other is out of [—2, 2], For the former, the corresponding eigenstates 
are propagating ones. For the latter, the wave functions are localized, which represent 
bound states. We have found that in our model of two coupled rings, there could exist 
bound states, which are localized around the junction. 

The potential applications of these localized states are unclear at this moment. We 
expect that if we couple many rings together and form a kind of nano-ring network, then 
at each junction there exists a localized state. The wave functions of the localized states 
at different junctions could have small overlaps, and so we expect that electrons are able 
to tunnel from one junction to a nearby junction. In this way, we have obtained a new 
lattice, formed by junctions between nearby nano-rings. Obviously, the lattice constant 
and the hopping constant can be adjusted by changing the size of the nano-rings. 

Finally, we would like to point out that the method we have used in this paper is not 
restricted to the single band model. It is also possible to be applied to multi-band tight 
binding models. In such circumstances, the transfer matrices are more complicated. 
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Appendix A. The calculation of L n 

In this appendix we derive the expression of L n , the result of which is shown in the main 
text from equation (|8|) to (12). 


From (j6j the transfer matrix L is 

A -1 
1 0 


L = 


(A.l) 


where A is a real number, corresponding to the eigenvalue of the original Hamiltonian 
h. 
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Let’s try to find the eigenvalue of L first. Assume (a,fc) T is an eigenvector of L 
with eigenvalue 77 . So 


A -1 
1 0 


b l =V \ b 


or 


= 0 


\ — Tj —1 

1 —7/ 1 \ b 

In order to get a non-zero eigenvector, we have to require 
A — 77 —1 


or 


Thus 


1 =0 

1 —77 

(A-77X-77) + 1 = 0 

A ± a/A 2 — 4 

V = -o- 


(A. 2 ) 

(A.3) 

(A.4) 

(A.5) 

(A. 6 ) 


Now according to the value of A, it is necessary to separate the discussion into three 
different situations 

i) if -2 < A < 2 

In this case we can set A = 2 cos#, where 0 < 9 < n. And so 


A ± a/A 2 - 4 ±i0 

Vl,2 = -o- = e 


For 771 = e , from (|A.2|) the corresponding eigenvector is 

<i\ ' ' ~ ie 

b\ 


(A.7) 


(A. 8 ) 


For 772 = e 1 , from (A.2) the corresponding eigenvector is 


02 

b 2 


1 

AO 


(A.9) 


Now we can do a liner transform to diagonalize L. The transformation matrix is formed 
as 

' 3 ie 1 


A = 


CL 1 Cl 2 


bi b 2 I \ 1 e 

and L can be written in the form 

L = ADA- 1 = A [ m 


id 


A 


-1 


(A. 10) 


(A. 11) 
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where 


D = 


Vi 


AO 


V-2 


-,—iO 


and 


A -1 = 


-Cl 2 


AO 


det(A) y —b\ a\ 

is the inverse of A. 

Now L n can be calculated easily 


e 2ie — 1 \ —1 e 


i o 


w /I —i 


L n = (AD A~ l )(ADA~ l ) ■■■(A DA~ l ) = AD n A 


Since 


We have 


D n = 


L n = 


Vi 


And 


V2 


e id 1 


i nO 


AnO 


AO 


3 2i0 - 1 1 


AO 


i nO 


-1 e 


i 0 


Do this matrix multiplication, and after simplification we obtain 

sin(n+l)6> sin nO 


L n = 


sin 6 sin 6 

sin nO _ sin(n-l)fl 
sin 6 sin 6 


(A.12) 

(A.13) 

(A.14) 

(A.15) 

(A.16) 

(A.17) 


ii) if A = ±2 


We only discuss the situation of A = 2 here. In this case 
A ± a/A 2 — 4 


n = 


= 1 


(A. 18) 


There is only one eigenvalue. From (A.2) we have the corresponding eigenvector 


a 

b 


(A. 19) 


There is no second linearly independent eigenvector. Thus L cannot be diagonalized. 
As a compromise, we try to find a vector (c, d) T that fulfils 



(A.20) 


(A.21) 


(A.22) 













19 


Now we form a transformation matrix 

a c \ l 11 


A ' b d 
and write L in the form 

L = ADA- 1 = A 


where 


and 


D = 


r] 1 
0 


1 0 


r] 1 

0 7] 

1 1 

0 1 


A 


-i 


A -1 — 


0 1 


1 -1 

is the inverse of A. 

Similar to the case —2 < A < 2, we have 
L n = AD n A ~ 1 

Now 


So 


D n = 


L n = 


rj 1 \ I r] n nr] n 1 


0 7] 

1 1 
1 0 


0 rf 


1 n 
0 1 


1 n 
0 1 


0 1 
1 -1 


71+1 

n 


— Tfl 

— (n-1) 


iii) if A > 2 or A < —2 

We only discuss the situation of A > 2 here. 
Set A = 2 cosh k, where k > 0. So 


+ 1,2 = 


A ± \/A 2 — 4 


= e 


±k 


For r/i = e k , from (A.2) the corresponding eigenvector is 

~./c 


a i 

b\ 


1 


For ij 2 = e k , from (|A.2|) the corresponding eigenvector is 


a 2 

b 2 


(A.23) 


(A.24) 


(A.25) 


(A.26) 


(A.27) 


(A.28) 


(A.29) 


(A.30) 


(A.31) 


(A.32) 
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Now we use a liner transform to diagonalize L. Form the transformation matrix 


A 


CL\ &2 

bi b 2 


1 


1 


and L can be written in the form 


where 


and 


L = ADA- 1 = A ^ >]1 ^ j A -1 




is the inverse of A. 


Then 


L n = AD n A ~ 1 


and we have 


So 


or 



Q—nk 



L n 


sinh(n+l)/c 
sinh k 
sinh nk 
sinh k 


_ sinh nk 

sinh k 

sinh(n— l)k 
sinh k 


(A.33) 

(A.34) 

(A.35) 

(A.36) 

(A.37) 

(A.38) 

(A.39) 

(A.40) 
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